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Abstract. The area distance to a convex plane curve is an important concept in computer vision. In this paper 
we describe a strong link between area distances and improper affine spheres. This link makes possible a better 
understanding of both theories. The concepts of the theory of affine spheres lead to a new definition of an area 
distance on the outer part of a convex plane arc. Also, based on the theory of discrete affine spheres, we propose 
fast algorithms to compute the area distances. On the other hand, area distances provide a good geometrical 
understanding of improper affine spheres. 
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1 Introduction 

The area distance of convex plane curves is an important concept in computer vision. These distances can be useful in 
matching two images of the same object obtained from different points of view [7]. It can also be seen as an erosion, a basic 
concept of mathematical morphology ([6]). 

Improper affine spheres are surfaces in whose affine normals at all points are parallel. In this paper, we point out the 
strong connection between these area distances and improper affine spheres. This connection is used in the development of 
the theory of area distances: based on the theory of affine spheres, we propose a new definition of area distance on the outer 
part of a convex curve and new algorithms for computing these distances. This link is also interesting from the point of view 
of the theory of improper affine spheres, since it provides a geometrical interpretation of them. 

Let us make this connection more precise: we first remark that the area distance / to a convex plane curve satisfies the 
Monge-Ampere differential equation det{D^{f)) — —1 with boundary conditions f — and V(/) = (see [8]). We show 
in this paper that the graph of / is an indefinite improper affine sphere, with strictly positive Pick invariant. On the other 
hand, we show also that, at least locally, any indefinite improper affine sphere is an area distance. 

For the outer part of the convex curve, we propose to define an area based distance by the Monge-Ampere differential 
equation det{D^{f)) = + 1 with boundary conditions / OandV(/) — 0. Then the graph of / is a definite improper affine 
sphere, also with strictly positive Pick invariant. We consider in this paper the case of an initial analytic curve. In this case, 
it is possible to find explicitly a solution to the Monge-Ampere equation, and also to describe its geometrical properties. 

Consider now a convex polygon as a discretization of the convex curve. By following the asymptotic lines, we propose a 
fast evolution algorithm that computes exactly the area distance. We also show that this exact area distance defines a discrete 
indefinite improper affine sphere, as defined in [5]. We remark that the method proposed in this paper differs completely 
from that of [9], since the latter considers curves defined in implicit form. For the outer part of the polygon, we also propose 
a fast evolution algorithm that computes a new distance. This graph of this new distance has the remarkable property of 
being a discrete definite improper affine sphere, as defined in [5]. 

In this context, there is a natural duality between points in the inner and the outer part of the convex curve. An interesting 
fact is that, although this duality is not area preserving, it preserves the measure J^^^dxdy, where J denotes the Pick invariant 
of the corresponding graph. From the discrete point of view, it is interesting to observe that a mesh with planar crosses in the 
inner part of the curve changes smoothly along the curve to a mesh with planar quadrilaterals in the outer part. 

This paper is organized as follows: In section 2, we review the basic concepts related to improper affine spheres, both 
smooth and discrete. In section 3, we review the definition of area distance in the inner part of a curve and show its strong 
link with indefinite improper affine spheres. In section 4, we propose the new definition of area distance in the outer part of 
the curve, show that its graph is an definite improper affine sphere and describe the duality between the inner and the outer 
area distances. 

Notation. For three vectors X, Y and Z in the space, denote by [X, Y, Z] the determinant of the 3x3 matrix whose columns 
are the vectors X, Y and Z. For two vectors X and Y in the plane, denote by [X, Y] the determinant of the 2x2 matrix 
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whose columns are the vectors X and Y . Also, denote by X* the transpose of the matrix X and by R the ninety degrees 
rotation in the anti-clockwise direction. Observe that [X, Y] = —X* ■ R-Y. 

2 Improper afflne spheres 

2.1 Basics 

For a surface S parameterized by q{u, v) = {x{u, v), y{u, v),z{u, v)) let L = qv,quu], M = [qu, Qv, Quv] and N = 
[Qu, Qv, Qvv]- We say that S is non-degenerate if LN — ^ 0. For a non-degenerate surface, the Blaschke metric is given 
by 

_ Ldu'^ + 2Mdudv + Ndv^ 

^~ |i7V-M2|l/4 • 

It is definite or indefinite according to LN — being positive or negative. In the definite case, the affine normal is defined 
as ^ = ^^{q), while in the indefinite case £, = — iA(g), where A((j') denotes the laplacian of each coordinate with respect 
to the Blaschke metric. 

An improper affine sphere S* is a surface whose affine normals at all points are parallel. We shall assume that the affine 
normals are parallel to the z-axis. Under this hypothesis S is locally the graph of a function z = f{x, y). The following 
proposition is well-known (see [1]). 

Proposition 1 S is an indefinite improper affine sphere if and only if det{D^{f)) — —c, and a definite improper affine 
sphere if and only ifdet{D^{f)) = c, where c is a positive constant. In both cases, the affine normal is constant and equal 
to (0, 0, c). 

In the case S is the graph of /, we shall write v) = {p{u, v),f{p{u, v))). The coefficients of the Blaschke metric 
can be calculated by the following lenoma, whose proof is a straightforward calculation: 

Lemma 2 

L = \pu,Pv]D'^{f){Pu,Pu) 
N= \pu,Pv]D'^{f){Pv,Pv) 

2.2 Asymptotic and isothermal directions 

For an indefinite Blaschke metric, one can find parameters (u, v) such that L = N = 0. These parameters are called 
asymptotic parameters and the corresponding tangent vectors are called asymptotic directions. Lines whose tangent vectors 
are asymptotic directions are called asymptotic lines. Lemma 2 shows that the projections of the asymptotic directions in the 
{x, y)-plane vanish the quadratic form D'^{f), and we shall also call them asymptotic directions. 
Using asymptotic parameters, the structure equations become 

Quu — ~^Qu "I" "^Qv 

quv = ~uj^ 

where Ot, = 5„ = and uj = [qu,qv,^]. By a good choice of the asymptotic parameters, we can make a and b constants. 
The Pick invariant is given by j = ^ (see [4] and [5]). The equations for the planar component are 

{Puu = ~^P'^ ^ TUP^ 
Pvv = "I" 'iJ'Pv 

Puv = 0, 

withw = = —fuv From these equations one obtains [pa-.puu] = o and = —b. So, if the Pick invariant 

does not vanish, the planar asymptotic lines are convex. 

For a definite Blaschke metric, one can consider also asymptotic parameters, but they are complex ([4]). Consider 
complex parameters z = s + it and z = s — it such that L{z, z) = N{z, z) = and iM {z,z) = CP. Then, in terms of s 
and t, L{s, t) = N{s, t) = iCP and M{s, t) =0 ([2]). We call such coordinates isothermal. 
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In order to differentiate from the indefinite case, we shall use capital letters P, Q, F, fl, A, B and J to describe the 
structure equations, as follows: 



where ft = —i[Qz, Qy, £,], A = i[Qz, Qzz,£,] and B = — Qyi, £,]■ By a good choice of the isothermal parameters, 
we can make A and B constants. The Pick invariant is given by J = ^ (see [2] and [5]). The equations for the planar 
components are 



2.3 Discrete improper afflne spheres. 

In [2], definitions of discrete proper affine spheres are proposed, both in the indefinite and in the definite case. In [5], 
these definitions are generalized to discrete improper affine spheres, indefinite and definite. We describe now these latter 
definitions, with a sUght modification in the definite case. Denote by Z the set of integers. 

Definition 1 A map q : is a discrete indefinite improper affine sphere if it has the following properties: 

1. For any € Z"^, the points q{i, j),q{i + ^, j), q(i — 1, j), Q{i,j — + ^) '^''e co-planar 

2. There exists a direction^ in such that, for any € Z"^, thevectorq{i,j) + q{i + l,j + l) — q{i + l,j) — q{i,j + l) 
is parallel to ^. 

Definition 2 A map Q : Z^ ^ R^ is a discrete definite improper affine sphere if it has the following properties: 

1. For any G Z"^, the points Q{i, j), Q{i + 1, j), Q{i, j + 1), Q{i + + 1) are co-planar. 

2. There exists a direction ^ in R^ such that, for any S Z^, the vector Q{i,j + 1) + Q{i,j — 1) + Q{i — Ijj) + 
Q{i + 1, j) — 4:Q{i,j) is parallel to ^. 

3 Inner Area distances 

In this section we review some properties of area distances and show the connection between area distances and affine 
spheres. We show that the graph of an area distance is an indefinite improper affine sphere and that, at least locally, any 
indefinite improper affine sphere is the graph of an area distance of some convex plane curve. 

Moreover, we show that the area distance of polygons define discrete indefinite improper affine spheres and how this 
fact can be applied to construct a very fast algorithm for computing area distances. 

3.1 Definition and properties of tlie area distance 

Consider a smooth convex curve C in the plane without parallel tangent lines. C can have 2, 1 or endpoints. Denote by D 
the plane region whose boundary is C and the curve(s) obtained from C by a similarity of ratio ^ based at each endpoint of 
C (see figure 1). 

A chord is a line segment connecting 2 points of C. For a point p ^ D, consider chords I passing through p that, together 
with C, bound regions Di. Denote by l{p) the chord such that the area of the corresponding region is minimum. The 
area distance function f{p) is then defined as half of the area of -D(()3). Sometimes we shall call this function inner area 
distance, since in section 4 we shall define another area distance. 

Any p E D is the mean point of the extremities of the chord l(p). This important property was proved first in [6] (see 
also [7]). Another important property is that V(/) (p) is orthogonal to the chord l{p), with half of the length of it (see figure 
2). The third important property is that det{D'^{f){p)) — —1, for any p e D. This last property was first proved in [8]. The 
next lennma and proposition describe the 2 latter properties with more details. 





with O = -i[Pz,Pz\. Also A = i[Pz, Pzz] and B = -i[P-, Pyz\- 



3 



(a) Domain for a curve with 2 
endpoints. 



(b) Domain for a curve with 1 
endpoint. 



(c) Domain for a curve without 
endpoints. 



Figure 1 : Domains for the inner area distance function. 




Figure 2: Gradient of the area distance. 
Lemma 3 Denote by C{u{p)) andC{v{p)) the extremities of the chord l{p). Considering u and v as functions ofp, we have 

^V"J — [C'(u),C'(v)] ^\") [C'(u)fi'(v)Y 

Proof. Since 2p = C{u{p)) + C{v{p)), 21 = C"(M)V(zt)* + C'{v)\'{vY where / is the identity 2x2 matrix. Multiplying 
by C'{vYR, one obtains the first formula. The second one is analogous. 

Proposition 4 The area distance f of a convex arc C satisfies the following formulas: 

1. The gradient of f is given by V(/) = R{\{C{v) — C(u))). 

2. C'{u) andC'iv) satisfy the equations D^{f){C'{u), C'{u)) = 0, D'^{f){C'{v), C'{v)) = Oand D^{f){C'{u), C'{v)) = 

[C'{u),C'{v)]. 

3. det{D\f)) = -1. 
Proof. 

1. By Green's theorem, one can write 4/(p) = JJ[C(s) — p, C'{8)\ds. Using the rules for differentiating integrals, one 
obtains 

4V(/)(p) = [C{v) - p, C'{v)]V{v) - [Ciu) - p, C"(«)]V(«) + r RC'is)ds, 

J U 

AV{f){p) - R{ ^c'{u),C'{v)] ^ + [C'{v),C'{u)] ^ + " 
Since, for any w, one can write w = j^i~f^^T^ C'{u) + 

[c>"c>)]'^'(^)' conclude that 
V{f){p)=R{\{C{v)-C{u))). 
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2. Differentiating one obtains 

2D^{f) = R{C'{u)V{uy - C'{v)V{vy), 



and so 

^'^•^^ = [C'{u),C'{v)] {RC'{u)C'{vrR + RC'{v)C'{uyR) . 

One concludes that 

C'{uYD\f)C'iu) = C'{vYD'{f)C'{v) = 

and 

C'{urD\f)C\v) = C'{vfD'{f)C'{u) = [C'{u),C'{v)]. 

3. Since we have assumed that there are no parallel tangents, property 2 guarantees that D'^(f ) is non-degenerate. Now 
C"(u)*L>2(/)C"(u) = implies that D^{f)C'{u) = XRC\u), for some A ^ 0. And since C {vf {f)C' {u) = 
[C'{u),C'{v)l\ = l. So D'^{f)C'{u) = RC'{u) and similarly D'^{f)C'{v) = -RC'{v). Hence [D'^{f)C'{u), D'^{f)C'{v)] = 
[RC'{u), -RC'{v)] = -[C'{u),C'{v)], which implies that det(£'2(/)) = _i. 



3.2 Area distances as improper affine spheres 

3.2.1 Consequences of proposition 4 

Proposition 4 impUes that the graph of the area distance / to a smooth convex arc C is an indefinite improper affine sphere. 
Also, the parameterization 

qiu,v) = {^iC{u)+Civ))Ji^iCiu) + C{vm. 

is asymptotic. Hence C'{u) and C'{v) are asymptotic directions and the asymptotic fines are obtained from C by similarities 
of ratio 1/2. Direct calculations show that v) = \ [C'{v), C'{u)]. Finally, f{p) represents the area of the region bounded 
by the two asymptotic fines that start at p and the curve C itself (see figure 3). 

If the parameterization of the curve C is by affine arc length, then a = b = —j. In this case the Pick invariant is 
j = 4:[C' (v) , C {u)]~^ . Thus it vanishes only at points p such that at least one of the endpoints of minimal chord l{p) 
belongs to a fine segment of the original curve C. In particular, if the curve is strictiy convex, the Pick invariant never 
vanishes. 

3.2.2 Some explicit formulas for calculating /. 

By the first item of proposition 4, formulas = fxXu + fyUu and fy = fxXy + fyUy imply that 



X (u) — X (v) x' {u) 
y {u)-y {v) y' (u) 

X (u) — X {v) x' {v) 
y(u)-y {v) y' {v) 

It is also interesting to consider coordinates (s, t) defined by 2s = u + v and 2t = u — v. In these coordinates, since 
9s= 9u + 9v and gt = Qu - Qv, we have 



fu {u, v) = - 



fv {u, v) = 



1 


x{s + t) 


— X {s 


-t) 


x' (s + + x' (s 


-t) 


4 


y{s + t) 


-y{s 


-t) 


2/'(s + t)+2/' (s 


-t) 


1 


x{s + t) 


— X {s 


-t) 


x' (s + t) - x' {s 


~t) 


4 


y{s + t) 


-y{s 


-t) 


y' (s + 1) - y' {s 


-t) 
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3.3 Examples 

In this subsection we give explicit examples of area based distance function of convex smooth curves. 

Example 1 Consider the parabola parameterized by C {r) = (r,^^- Integrating [f^^ fv) = \{{u — v)"^ ,—{u — vY) one 
obtains f{u, v) = ^{u — v)^. Thus a parameterization of the affine sphere in asymptotic coordinates is given by 

U + V u'^ +v'^ 1 3 

m,v) = {^-, — - — ,^{u-v) ) 

with u > V. Andiv{u, v) = j{u — v). In (s, t) coordinates. 



g2 _|_ ^2 ^3 

Q{s,t) = {s, — - — ,-) 



and uj{s, t) = ^t. 



Example 2 Consider the circle parameterized by C{r) = (cos(r), sin(r)). Although this curve admits parallel tangents, the 
above calculations work well, except at the center of the circle. Integrating (/„, /„) = i(— cos (m — w) + 1, cos {u — v) — l) 
one obtains f{u, v) = j{u — v — sm{u — v)). Thus a parameterization in asymptotic coordinates is given by 

q{u, v) = (i(cos(M) + cos(t')), i(sin(u) + sin(ti)), ^{u ~ v — sm{u — v))) 

with u> V. AndLj{u, v) = j sin(u — v). In (s, t) coordinates. 



q{s, t) = (cos(s) cos(f), sin(s) cos(t), i(2f — sin(2t))) 



andLo{s,t) = ^sin{2t). 



Example3 Consider the hyperbola parameterized by C{r) = (exp(r), cxp(— r)). Integrating {fu,fv) = j{—'2 + exp{u — 
v) + exp(w — m), 2 — exp(M — v) — exp{v — u)) one obtains f{u, v) = ^{v — u) + j{exp{u — v)— exp(u — u)). Thus the 
asymptotic parameterization is 

q{u, v) = -(exp('u) + exp(t;),exp(— u) + exp(— t;), (v — u) + 2 (^■''P (u — v)— exp {v — u))). 

with u > V. AndLv{u, v) = j{exp{u — v) — exp{v — u)). In (s, t) coordinates, 

q{s,t) = ^{exp{s + t)+exp{s-t),exp{-{s + t))+exp{t - s),-2t+^{exp{2t) -exp(-2t))). 

anduj{s,t) = |(exp(2t) — exp(— 2t)). 

Example 4 Consider the cubic parameterized by C{r) = (r, r^). Although this is not a convex curve, the above cal- 
culations can be done. Integrating the vector field {fmfv) = ji'^u^ ~ Su^ii + v'^,u'^ — one obtains 
f(u, v) = \{u — v)^ (u + v). Thus the asymptotic parameterization is 

q{u,v) = ^((u + v),u^ + v^, ^ (w - v)^ {u + v)). 



with u> V. Anduj{u, v) ~ — v^). In (s, t) coordinates, 

q{s,t) = {s,s{^t^ + s^),2t^s) 

andu{s,t) = 3st. 
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Example 5 Consider the curve parameterized by C{r) = (r, r^). Integrating the vector field (/„, fy) = j{3u'^ — 4:U^v + 

l/3ul _ 3t£ 
4^ 5 5 



t;*, Av^u — — M^) one obtains f{u, v) = u^v + uv'^). Thus the asymptotic parameterization is 



Qiu, V) = -{u + v,u +v , -{— u v + uv )). 

with u> V. And u{u, v) = — v^. In (s, t) coordinates, 



q{s, t) = {s, t^ + s^ + 6sH'^,4t^s'^ + ^) 



andw{s,t) = 2t^ + 6ts^. 



3.4 Local characterization of indefinite improper affine spheres 

In this subsection, we show that locally, and up to a constant, any indefinite improper affine sphere with non-zero Pick 
invariant is the graph of the area distance of a smooth convex plane curve C. More precisely, we have the following theorem: 

Theorem 5 Let U be an open domain in the (?/, v)-plane whose closure U is contained in the domain of the asymptotic 
parameterization of an indefinite improper affine sphere S. Assume that, restricted to U, S is the graph of a function f{p). 
Then there exists a convex curve C in the plane and a constant K such that f + K is the area distance ofC. 

3.4.1 Some properties of indefinite improper affine spheres 

Consider an indefinite improper affine sphere with strictly positive Pick invariant j. Assume, w.l.o.g., that the affine normal 
is ^ = (0, 0, 1) and consider that S is the graph of a function /. 

Lemma 6 The following properties hold (see figure 3): 

1. We have that = Rpu and D'^{f){py) = -Rpv 

2. Let TTi{p) = p + i?V(/)(p) and tt2{p) = P — Then 7ri(p) is constant along an asymptotic line v = vq and 
772 (p) is constant along an asymptotic line u = Uq. 

Proof. 

1. Since is non-degenerateand Z)^(/)(p„,p„) = 0,D^(/)pu = Ai?p„, forsome A 7^ 0. Andsince£>^(/)(p„,p„) = 
u){u, v) = \pu,Pv]f A = 1. A similar reasoning shows that D'^{f){py) = —Rpy. 

2. Just observe that p„ + RD'^ (/)p„ = and p„ - RD"^ [f)p^=Q. 









-^m^ — V 










\^ f(p)=area 









Figure 3: The projections tti and 772. 
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3.4.2 Proof of theorem 5 

Let TTi and Ti2 be the projections defined in lemma 6. We have that Ci = T^iiU) and C2 — t^2{U) are compact arcs. It is not 
difficult to obtain a smooth convex arc C3 such that the concatenation C of Ci, C3 and C2 is smooth nd convex. Denote by 
g the area distance function associated to C. 

One can easily see that p = i(7ri(p) + 7r2(p)) and V{f){p) = \R{'K2{p) — ^i(p)). for any p G [/. So V(5)(p) = 
V{f){p), for any p . This implies that f — g is constant, which proves the theorem. 

3.5 Area distances to polygons 
3.5.1 Asymptotic grids 2d and 3d 

Let C be a convex polygon with vertices c;, < i < TV — 1. Denote its sides by the vectors 2Li — c^+i — and assume 
that [Li , Lj] > 0, for any j > i. Assuming j > i, define the following grid p{i,j) on the plane by 

/ ■ ■^ + Cj 
Pihj) = 

Note that 

+ = ^2L,^p{i + l,j + l)-p + 1) , 

so the grid is formed by parallelograms whose areas will be denoted by 

aij = b(« + - pihj) ,p{i-j + 1) -P{hi)] = [Li,Lj\ . 

(see figure 4). Define also 

i<k<l<j-l i<k<l<j~l 

Note that / (i, j) = if j = i or j = i + L 




p(4,4) p(4.5| p(5,5| 



(a) Points and parallelograms. (b) Areas defined by the parallel- 

ograms. 

Figure 4: Planar mesh of a discrete inner distance. 



Proposition 7 The map q (i, j) = [p {%,]) , / (i, j)) is a discrete indefinite improper affine sphere. 
Proof. 
L Let 

k — i 
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Then, since [Li, Li] = [Lj, Lj] = 0: 

Q{i,j + i) -QiiJ) 



Li-i, - ^ [Li-i, Li] 1 = (Z/j-i, [L, ij-i]) 
^ l=i ) 

L„- \ ={Li\LM) 

' i-i \ 
L„^[L,,L,] U(i,-,[i,i,-]) 



so whatever linear dependence is satisfied by Lj-i, and Lj will also be satisfied by the z coordinates of the three 
above vectors. This shows that g (z + 1, j) is in the same plane as g (i — 1, j), q — 1) and q (i, j). Similarly, one 
can show that q{i + l,j + 1) is also in this plane (see figure 5). 



2. From the equations above, note that 

q{i + l,j + l)-q{i,j + l) 

Qii + 'i-J) -Q{i,j) 

Subtracting, 



Li, 
Li, 



Lk, Li 

_k—i 

i-1 

Lki Li 



qii + l,j + l)-q + l)-qii + l,j) + q = (0, [L„ Li]) = (0, -aij) 
and hence these vectors are all parallel to the z-axis. 




Figure 5: The four line segments in bold are co-planar. 



3.5.2 Fast algorithm 

Now, if we define the level of a point to be fc = j — i we can actually obtain our grid in levels starting from the polygon. At 
level 0, p [i, i) = Ci and / (i, i) ~ 0. At level 1, p (z, z + 1) = i (q + q+i) and / (z, z + 1) = 0. At level (j - z) + 1, 

p{i,j + 1) =p{i,j)+p{i + l,j + l)-p{i + l,j) 
f + 1) = / + fii + l,j + l)-f{i + + [Li, Lj] 

This gives us a fast algorithm to calculate all the grid points. In each parallelogram, we can calculate the exact distance 
by a bilinear interpolation. In figure 6, one can see the result of this algorithm apphed to a polygon inscribed in an elhpse. 
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(a) A discrete improper afSne 
spliere. 



(b) Planar cross in bold. 



Figure 6: The inner area distance of a polygon inscribed in a circle. 



4 Outer area distances 

In this section, we consider the question of extending the area distance to a neighborhood EofCin the outer part of C. The 
idea is to solve the Monge- Ampere differential equation 

det{D^F){P) = +l,P€E 
VF(P) = 0,P G C 
F{P) = 0,PeC 

and define the area distance at P by F{P). It is clear that the graph of F defines an improper definite affine sphere. 

In this section, we shall first describe a solution to this problem in the case of an analytic curve. Then we prove some 
properties of this solution, including a relation with the area distance in D. Finally, we indicate how to obtain a discrete 
definite improper affine sphere as an outer area distance of a polygon. 



4.1 Definite improper affine splieres from analytic curves 

Assume that the parameterization C(r) = {x{r), y{r)) is analytic so that we can evaluate its coordinates for complex values 
of the parameter, z = s + it, and z = s — it. Since C is analytic and real on the real hne, the expression P{z,J) = 
^{C{z) +C{z)) is real. It represents the planar coordinates of the definite improper affine sphere defined below. 

Lemma 8 Consider the inner area function f in coordinates (s, t), and let F{s, t) = if{s, it). Then F{s, t) is real and 
V(F) = iR[£M}^Ml), 

Proof. From formulas (1) and (2), we can write 



Fs{s,t) 
Ft{s,t) 



x{z) — X (z) x' [z) + x' {z) 
y {z)-y (z) y' (z) + y' (z) 

X (z) — x (z) x' (z) — x' (z) 
y {z)-y (z) y' (z) - y' (z) 



(3) 
(4) 



We conclude that F{s, t) satisfies Fg = F^Xs + Fyys and Ft = F^xt + Fyyt, where F^ = -|(y(z) - yiz)) and Fy = 
- x{z)). Thus V(P) 



iR{C{z) — C{z)) is real. So F is also real. 



c(£)_c(£) -J Qi^sgj^vg aJsQ (j^at, since C{z) is analytic and real on the real line, V(F) = 



Proposition 9 The parameterization Q{s, i) = {P{s + it, s — it), F{s, t)) is isothermal and defines an improper definite 
affine sphere S. The surface does not depend on the choice of the analytic parameterization of the curve C. Also Q.{s, t) = 
—iio{s,it), A{s,t) = ia{s,it) and B{s,t) = —ib{s,—it). Hence the Pick invariant does not vanish for strictly convex 
curves C. 
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Proof. Since \7{F) = iR{^){C{z) — C{z)), the same proof of proposition 4, item (3), with {z,z) in place of (u, v) imphes 
that det{D'^{F)) = +1 and that Q{s, t) is an isothermal parameterization of a definite improper affine sphere. If we begin 
with a different analytic parameterization C{tp{z)), denote by F\ the corresponding function. From the above formulas, 
we have V(-Fi)(^;) = V {F){ij;{zy)ij:' {z), which imphes that F\{z) = F{ip{z)). Thus the surface is independent of the 
parameterization . 

SmceQ,{z,z) = {[C {z) , C {z)] = i[P^, P^-], we obtain f2(s,t) = ^[Pt,Ps]. On the other hand, w(s,t) = ^[p8,Pt]- 
Thus ri(s,t) — —iuj{s,it). 

Finally, ^(s, i) = ^[C'{s + it),C"{s + it)] =ia{s,it) and B{s,t) = -^[C'{s - it),C"{s - it)] = -ib{s,-it). 



4.2 Examples 

Given a curve C, assume that we can find a parameterization C(s) = {x{s),y{s)), with x{z) and y{z) analytic functions. 



~2 I —2 



Example 6 Consider the parabola C{r) = (r, ^) of example 1. Since X{z, z) = and Y{s, t) = ^ , we obtain 

2 



X{s,t) = sandY{s,t) = Also, F{s,t) = if{s,it) = ^. So 



- 1^ t^ 
Q{s,t) = {s,^—,-) 

t> Q, defines a definite improper affine sphere. The area element of the Blaschke metric is Q.{s, t) = —iu}{s, it) = |. 

Example 7 Consider the circle C{r) = (cos(r), sin(r)) of example 2. Since X(z,z) = '^°^(^)+™^(^) q^^^ Y{z,z) = 
sin(z)+sin(?) ^ ^^^^.^ ^^^^ ^ ^^^^^^ cosh(t) and Y{s, t) = sin(s) cosh(t). Also F{s, t) = if{s, it) = - t). 

Hence 

1 sinh(2t) 

Q{s,t) = (cos(s) cosh(t),sin(s) cosh(t), -( t)) 

t > 0, defines a definite improper affine sphere. The area element of the Blaschke metric is f2(s, t) = —iuj{s, it) = 

sinh(t) cosh(t) 
2 



Examples Consider the hyperbola C{r) — {cxp{r),cxp{—r)) of example 3. Since X {z, z) = '^^p(^)+'^^p(^) andY{z,z) 
Hence 



2 

we obtain X{s, t) = exp(s) cos(t) and Y{s, t) = exp(— s) cos(t). Also, F{s, t) = if{s, it) = t — SiSiEii 



Q{s, t) = (exp(s) cos(t), exp(— s) cos(t), t — £1^^^^) 
with t > 0, defines a definite improper affine sphere. The area element of the Blaschke metric is given by fl{s, t) = 

Example 9 Consider the cubic C{r) = (r, r^) of example 4. Since X{z,z) = ^{z + z) and Y{z,z) = ^{z^ + z^), we 
obtain X{s, t) = s, Y{s, t) = - Sst"^. Also F{s, t) = if{s, it) = 2t^s. Hence 

Q{s,t) = {s,s^ -3st'^,2t^s) 

t > 0, defines a definite improper affine sphere. The area element of the Blaschke metric is Q.{s, t) = —iuj{s, it) — 3st. 

Example 10 Consider the curve C{r) = {r,r^) of example 5. Since X{z,z) = and Y{z,z) = ^ we obtain 
X{s, t) = s and Y{s, t) = s'^ - 6sH'^ + t^. Also F{s, t) = if{s, it) = 4sH^ - ^t^. Hence 

Q{s, t) = {s, - 6sH^ + ^^ AsH^ - ^i^) 

5 

t > 0, defines a definite improper affine sphere. The area element of the Blaschke metric is Q{s,t) = —ico{s, t) = 6s'^t — 2t^. 
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4.3 Isothermal tangent lines to the curve 

4.3.1 Another change of variables 

Consider new coordinates {u, v) defined by s + t = u and s — t = v. 

Lemma 10 Fix a point P{uq, vq), uq > vq and consider the lines P{u, uq), < u < uq and P{uo, v), vq < v < uq. 
These lines touches tangentially the curve at the points P{vo, Vq) and P{uo, Uq), respectively. 

Proof. Straightforward calculations shows that 

Pu = l{il + i)C'{z) + {l-i)C'{z)) 

and 

P. = \{{l-i)C'{z) + {l + i)C'{z)). 

Thus, at P{uo, Wo), Pv = Similarly, at P{vq, vq), P„ = 

We shall call the above lines the isothermal tangent lines. It is interesting to observe that in the case of a parabola, these 
isothermal tangent lines are in fact straight lines (see example 6). Also, in examples 6, 7, 8 and 9, the isothermal lines starting 
at a point P{uo, vq) do not meet C before the tangency points P{uo, uq) and P{vo, vq). We shall refer to this property as 
non-crossing isothermal tangents. Example 10 does not have this property, i.e., the isothermal tangent lines meet C before 
they arrive at the tangency points. 

4.3.2 Geometric interpretation of F as an outer area distance 

We shall assume from now on that the tangent isothermals are non-crossing. Considering coordinates {s,t), we can define 
a bijection between the iimer and the outer parts of C, by corresponding the points p = p{s, t) and P = P{s, ti). This 
correspondence can be defined more geometrically as follows: For p G D, consider the asymptotic lines that passes through 
p and denote by C (u) and C (v) the points where they touch tangentially C. Consider then the tangent isothermal Unes in E 
that touch tangentially C at C{u) and C{v). The intersection of these fines is P (see figure 7). 



C(v) 




C(u) 



Figure 7: The correspondence between p and P. 

For a point P = P{uo, vq), denote by A{P) the region bounded by the isothermal tangent lines and the part of the curve 
C(r) with vo < r < uo (see figure). Next lemma shows that the area of A{P) is F{P)/2. This property justify the name 
outer area distance for the function F. 

Lemma 11 Assume that the isothermal tangent lines starting at P = P{uo, vq) are non-crossing. Then the area ofA{P) is 
F{P)/2. 

Proof. Straightforward calculations shows that [P„, P-^] = ^. Thus the area of A{P) is given by half of the integral of O 
over the triangle T{uo, vq) whose vertices are (mq, vq), {uq, uq) and {vq, vq). The corresponding inner region a{p) has area 
equal to the integral of w over the same triangle. And we know that this area is equal to f{p). Thus 

F{P)=if{p) = i / co{u,v)dudv = / rt{u,v)dudv 

Jt{uo,vo) JT{uo,va) 

and so the area of A{P) is F{P)/2. 
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4.3.3 Relation between the indefinite and definite affine splieres 
Proposition 12 Let b C D and B c E be corresponding regions. Then 

Jq(b) JQ{B) 

where the integrals are taken with respect to the Berwald-Blaschke metric. 

Proof. Remember that j = abw~^ and J = ABQ,~^. We can assume w.l.o.g. that the curve C is parameterized by affine 
arc-length. Then j = J = 1/16. Up to these constants, the above integrals correspond to the areas, in the (u, w)-plane, of 
the regions that represent b and B, respectively. Since, in {u, v) coordinates, both regions are the same, the proposition is 
proved. 

4.4 Discrete outer area distances to polygons 

Assume that C is a polygon with vertices Pi = {xi,yi), 1 < i < N. In order to mimic the continuous case, we must consider 
that X and y are functions of a variable u € Z and then extend these functions to discrete complex analytic functions. 

There are several definitions of discrete analytic functions. We shall adopt a classical one (see [3], ch.5). Consider g : 
Z"^ iJand h : (Z^)* — » R, where (Z^)* denotes the dual lattice, g and h are complex conjugates if g{u+l, v)—g{u, v) = 
h{u + 1 /2,v + I /2) - h{u+ I /2,v- 1/2) saAg{u,v + l)- g{u,v) = -{h{u + l/2,v+l/2) -h{u-l/2,v+l/2)),for 
any (w, v) E Z"^. These equations are called discrete Cauchy-Riemann equations, and we say that g + ihis discrete analytic. 

Using the above definition, we can extend x and y to discrete analytic functions X + ih\ and Y + ih2. In fact, these 
functions are uniquely defined if we consider that hj{u + 1/2, 1/2) = —hj{u + 1/2, —1/2), j = 1, 2. This condition is 
natural for any analytic function that is real on the real line. 

As in the continuous case, we must find F(u, v), {u, v) E Z^ such that (/12, —hi) = V(F). We write 

F{u+l,v)- F{u, v) = h2{u + 1/2, V + 1/2) ■{X{u + 1, v) - X{u, v)) -hi{u + 1/2, v + 1/2) • {Y{u + l,v)- Y{u, v)) 
and 

F{u, v + 1)- F{u, v) = h2{u + 1/2, t; + 1/2) • {X{u, v + 1)- X{u, v)) -hi{u + 1/2, t; + 1/2) • {Y{u, v + 1)- Y{u, v)) 
and, using the discrete Cauchy-Riemann equations, we obtain 



F{u+l,v) - F{u,v) = 
F{u,v + 1)- F{u,v) = 



hi{u+l/2,v-l/2) hi{u+l/2,v + l/2) 
/i2(m+ 1/2,'y- 1/2) /i2(m + 1/2,1; + 1/2) 

hi{u-l/2,v + l/2) hi{u + l/2,v+l/2) 
h2{u-l/2,v + l/2) h2{u + l/2,v+l/2) 



(5) 
(6) 



In order to simplify notations, we shall denote by h the vector {hi ,h2)- 
Lemma 13 There exists F : Z^ ^ R satisfying equations (5) and (6). 

Proof. One has to show that the discrete derivative DFu(v) of the right hand of (5) with respect to v is equal to the discrete 
derivative DFv{u) of the right hand of (6) with respect to u. We have 



DFu{v) = [ h{u + l/2,v- 1/2) + h{u + l/2,v + 3/2) h{u -h 1/2, u -h 1/2) ] 



and 



DFv{u) = [ h{u + 1/2, v + 1/2) h{u - 1/2, u -I- 1/2) + h{u + 3/2, u + 1/2) ] 
Thus DFv{u) — DFu(v) is given by 

[ h{u + l/2,v + 1/2) h{u - 1/2, u -h 1/2) -h h{u + 3/2, v + l/2) + h{u + l/2,v- 1/2) -|- h{u +l/2,v + 3/2) ] 

which is equal to zero, since hi and /12 are discrete harmonic. 
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We shall denote P = {X,Y) and Q = {X,Y,F). Also, define the co-normal vector hy u = (— /i2,/ii,l) = 
(-V(i^), 1). We can write 

Q{u+l,v) -Q{u,v) = v{u+l/2,v+l/2)xv{u + l/2,v-l/2) (7) 
Q{u,v + 1)-Q{u,v) = -iy{u + l/2,v+ 1/2) xiy{u-l/2,v+ 1/2). (8) 

It follows directly from these equations that the quadrangles whose vertices are Q{u, v), Q{u + 1, v), Q{u, v + 1) and 
Q{u + 1, w + 1) are planar. In fact, each edge is orthogonal to v{u + l/2,ti + l/2). Thus Q{u, v) defines a conjugate net (see 
definition in [3]). Denote by /^{g) = g{u +l,v)+ g{u, v + l)+ g{u —l,v)+ g{u, w — 1) — 4:g{u, v) the discrete laplacian 
of g. It is clear that A{Q) is parallel to the z-axis. This follows directly from the fact that X and Y are discrete harmonic, 
i.e., A(X) = A{Y) = 0. Thus Q{u, v) defines a discrete definite improper affine sphere. 

We can also calculate A(F). Denote the vectors P{u + 1, w) - P{u, u), P{u, v + 1) — P{u, u), P{u - l,v) - P{u, v) 
and u — 1) — v) by Vi,V2, and v^, respectively. Since v) is discrete harmonic, vi + f 2 + I's + ^'4 = 0. 
And the area A{u, v) of the quadrangle whose vertices are P{u + l,v), P{u, v + 1), P{u — 1, v) and P{u, v — 1) is given 

by [V3,V2] + [VI,V4]. 

Lemma 14 

A{F){u,v) = A{u,v) 

Proof. We take as reference the vector h = h{ii, + 1/2, f + 1 /2). Then 

h{u-l/2,v + l/2) = h + V2 
h{u + 1/2, V -1/2) = h-vi 
h{u-l/2,v-l/2) = h + V2 + V3 



So 



= h — vi — V4 



Fiu+l,v) - F{u,v) = -[h~vi,h] = [vi,h] 

F{u,v + 1)- F{u,v) = [h + V2,h] = [v2,h] 

F{U-1,V) - F{U,V) = [h + V2+V3,h + V2] = [v3,h] + [v3,V2] 

F{u,v — 1) — F{u,v) = —[h — vi—V4,h — vi] = [v4,h] + [vi,V4] 



We conclude that 



F{u + l,v) + F{u -l,v)- 2F{u,v) = [v3,V2] + [vi,V4\ = A{u,v), 
thus proving the lemma. 




Figure 8: Outer area distance of a polygon inscribed in a circle: observe the planar quadrilateral in bold. 



One can see in figure 8 the outer area distance to a polygon inscribed in a circle. In our discrete construction, it is 
worthwhile to observe an asymptotic net and a conjugate net meeting along a curve. Unfortunately we were not able to give 
a geometric interpretation of the function F in the discrete case as we have done in the continuous case. 
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5 Conclusions 

We have shown a very close connection between area based distance, a widely used concept in computer vision, and improper 
affine spheres. Both theories may take benefit from this connection. 

From the point of view of the theory of area based distances, this link allow us to propose a new area based distance 
outside a convex region and to develop fast algorithms for computing the inner areas. 

From the point of view of the theory of improper affine spheres, the approach give a very geometrical description of 
these surfaces, both in the smooth and in the discrete case. 
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(a) Polygon inscribed in a cir- 
cle: (u, v) coordinates inside and 
(s, t) outside. 



(b) A closer look. 




(c) Polygon inscribed in a (d) A closer look, 

circle: (m, 1)) coordinates inside 
and {u, v) outside. 




(e) Polygon inscribed in a (f) A closer look, 

parabola. 

Figure 9: Internal asymptotic and external isothermal lines of polygons. 
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